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Abstract 


In this paper, we consider a coupled PDE system describing phase separation and dam¬ 
age phenomena in elastically stressed alloys in the presence of inertial effects. The material 
is considered on a bounded Lipschitz domain with mixed boundary conditions for the dis¬ 
placement variable. The main aim of this work is to establish existence of weak solutions 
for the introduced hyperbolic-parabolic system. To this end, we first adopt the notion of 
weak solutions introduced in |HK11 |. Then we prove existence of weak solutions by means 
of regularization, time-discretization and different variational techniques. 
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1 Introduction 

In micro-electronic materials such as solder alloys, different physical processes are shaping the 
micro-structure. For a realistic description of these structures, phase separation, coarsening and 
elasticity as well as damage phenomena have to be taken into account. A fully coupled system 
has been originally studied in |HK11) and further developed in [HK13b] allowing, for instance, 
inhomogeneous elastic energy densities. The corresponding degenerating case has been analyzed 
in [HK12] . To the authors’ best knowledge, before these works, phase separation and damage 
processes have only been investigated independently of each other in the mathematical literature. 

Phase separation and coarsening phenomena are usually described by phase-held models of 
Cahn-Hilliard type. The evolution is modeled by a parabolic diffusion equation for the phase 
fractions. To include elastic effects, resulting from stresses caused by different elastic proper¬ 
ties of the phases, Cahn-Hilliard systems are coupled with an elliptic equation in the case of a 
quasi-static balance of forces. Such coupled Cahn-Hilliard systems with elasticity are also called 
Cahn-Larche systems. Since in general the mobility, stiffness and surface tension coefficients 
depend on the phases (see for instance [BDM07| and [BDDM07] for the explicit structure de¬ 
duced by the embedded atom method), the mathematical analysis of the coupled problem is very 
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complex. Existence results were derived for special cases in |CMP00llGar00[rBP05) (constant mo¬ 
bility, stiffness and surface tension coefficients), in |BCD~*~0^ (concentration dependent mobility, 
two space dimensions), |SP13bl ISPlSa] (concentration dependent surface tension and nonlinear 
diffusion) and in |PZ08) in an abstract measure-valued setting (concentration dependent mobility 
and surface tension tensors). 

Damage behavior, however, originates from breaking atomic links in the material from a 
microscopic point of view whereas a macroscopic theory may specify damage in the isotropic case 
by a scalar-valued variable related to the proportion of damaged bonds in the micro-structure 
of the material with respect to the undamaged ones. According to the latter perspective, phase- 
field models are quite common to model smooth transitions between damaged and undamaged 
material states. Such phase-field models have been mainly investigated for incomplete damage 
which means that damaged material cannot loose all its elastic energy. 

Existence and uniqueness results for damage models of viscoelastic materials are proven in 
[BSS05] for scalar-valued displacements. Higher dimensional damage models are analytically in¬ 
vestigated in |BS041 [MR061 IMTToI IKR Zl^ IR R 12] and, there, existence and regularity properties 
are shown. A coupled system describing incomplete damage, linear elasticity and phase sepa¬ 
ration appeared in [HKlll IHK13b) . There, existence of weak solutions has been proven under 
mild assumptions, where, for instance, the stiffness tensor may be material-dependent and the 
chemical free energy may be of polynomial or logarithmic type. All these works are based on 
the gradient-of-damage model proposed by Fremond and Nedjar |FN96] (see also |Fre02] 'l which 
describes damage as a result from microscopic movements in the solid. The distinction between 
a balance law for the microscopic forces and constitutive relations of the material yield a satis¬ 
fying derivation of an evolution law for the damage propagation from the physical point of view. 
In particular, the gradient of the damage variable enters the resulting equation and serves as a 
regularization term for the mathematical analysis as well as it ensures the structural size effect. 
Internal constraints are ensured by the presence of non-smooth operators (subdifferential opera¬ 
tors) in the evolution system. Hence, in the case that the evolution of the damage is assumed to 
be uni-directional, i.e. the damage process is irreversible, the microforce balance law becomes a 
doubly-nonlinear differential inclusion. 

The main aim of this paper is to generalize the results for hyperbolic-parabolic damage 
systems introduced in [HK13a] to coupled phase-field systems describing phase separation and 
damage processes in the presence of inertial terms with mixed boundary conditions on non¬ 
smooth (Lipschitz) domains. The novelty of this contribution is to obtain existence results for 
phase separation with elasticity including inertial effects and damage processes on Lipschitz 
domains. We first utilize and adjust the notion of weak solutions introduced in |HK11| . Then, 
we prove existence of weak solutions by means of regularization, time-discretization and different 
variational techniques. To this end, an energy estimate has, for instance, to be established and 
several convergence properties are shown. 

1.1 Energies and evolutionary equations 

Here, we qualify our model formally and postpone a rigorous treatment to Section The 
presented model is based on two functionals, i.e. a generalized Ginzburg-Landau free energy 
functional £ and a damage pseudo-dissipation potential TZ (in the sense by Moreau). The free 
energy density g} of the system is given by 

V?(e(u),c, Vc,z, V 2 ;) := -\Vz\p + ^\Vc\'^ -\-W{c,e{u), z) f{z) 'S{c), (1) 

p 2 

where the gradient terms penalize spatial changes of the variables c and z. W denotes the 
elastically stored energy density accounting for elastic deformations and damage effects, / is the 
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damage dependent potential and stands for the chemical energy density. 

The overall free energy £ of Ginzburg-Landau type has the following structure: 


£{u,c,z) := z,\'z) + I[o^oo){z)j dx. (2) 

In this context, d[o.oo) signifies the indicator function of the subset [0,oo) C R, i.e. /[o,oo)(a:) = 0 
for X £ [0,oo) and /[o_oo)(2:) = oo for a; < 0. We assume that the energy dissipation for the 
damage process is triggered by a rate-dependent dissipation potential TZ of the form 


TZ{z) 




0 



dx. 


(3) 


The governing evolutionary equations for a system state q = {u, c, z) can be expressed by 
virtue of the functionals @ and More precisely, the evolution is driven by the following 
hyperbolic-parabolic system of differential equations and differential inclusions: 


diffusion: 

Ct = div(m(c, z)V/r), 

(4a) 


P = -Ac -1- W^c{c, e{u), z) + 4''(c), 

(4b) 

balance of forces: 

utt - div (We(c, e{u),z)) = 1, 

(4c) 

damage evolution: 

0 £ dz£(u, c, z) -I- dzTZ{dtz) or equivalently 

(4d) 


zt - ApZ + W^z{c, e(u),z) + f'{z) ^ = 0, 

(4e) 


f £ 9/[o,oo)(^)i 

(4f) 



(4g) 


The Cahn-Hilliard system (Ha1) - (l4bl) describes phase separation phenomena in alloys, the hyper¬ 
bolic equation formulates the balance of forces including inertial effects and the inclusion 
(l4dl)-(|4g| is an evolution law for the damage processes. The sub-gradients correspond to the con¬ 
straints that the damage is non-negative and irreversible. Let us note that linear contributions 
in / model damage activation thresholds. 

We choose Dirichlet conditions for the displacements u on a subset T of the boundary dft 
with > 0. Let & : [0,T] x T —>■ R" be a function which prescribes the displacements 

on r for a fixed chosen time interval [0,T]. The imposed boundary and initial conditions and 
constraints are as follows: 


boundary displacements : 
initial concentration : 
initial displacements : 
initial damage : 


M = & on Ld X (0,T), (5a) 

c(0) = c° in fl, (5b) 

m( 0) = 'Ut(O) = in fl, (5c) 

z(0) = in £L. (5d) 


Moreover, we use natural boundary conditions for the remaining variables on (parts of) the 
boundary: 


We(c, e(u), z) • 1 /= 0 onrNx(0,T), (6a) 

Vc-V = Vz ■ V = m{c,z)Vp ■ V = Q on (6b) 

where v stands for the outer unit normal to dVL. 

We like to mention that mass conservation of the system follows from the diffusion equation 
(Hal) and (l6bl) . i.e. 

f c{t) — c° dx = 0 for all t £ [0, T], 

Jn 
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In the next section, we state the precise assumptions that are needed for a rigorous analysis. 
Section [3] presents the main results. We give a notion of weak solutions evolved from |HK13a] 
and state the existence theorem in Subsection EH Since the proof is based on regularization 
techniques, we also give the weak notion and the associated existence result for the regularized 
system in Subsection EH In the main part, Section 01 the existence proof is carried out first for 
the regularized case and then for the limiting case. 

2 Notation and assumptions 

Throughout this work, let p > n be a constant and let C R" (n = 1,2,3) be a bounded 
Lipschitz domain. For the Dirichlet boundary Fd and the Neumann boundary Fn of we 
adopt the assumptions from [Berllj . i.e., Fd and Fn are non-empty and relatively open sets in 
dQ with finitely many path-connected components such that Fd D Fn = 0 and Fd U Fn = dfl. 

The considered time interval is denoted by [0, T] and fit := fl x [0, t] for t G [0, T], The partial 
derivative of a function h with respect to a variable s is abbreviated by h^s. The set {u > 0} 
for a function v G has to be read as {cc G fl | v{x) > 0} by employing the embedding 

^ C{Q) (because p > n). 

The elastic energy density W is assumed to be of the form 

Wic,e,z) = ^C{z){e-e*{c)) : (e-e*(c)), (7) 

where e* denotes the eigenstrain and C the material stiffness tensor which depends on the damage 
variable. For e*, we assume the linear relation e*(c) = ce with e G (Vegard’s law). We 

choose the stiffness tensor function C G C^([0,1]; £sym(R”^")), where £sym(R”^") denotes the 
linear mappings from into R”^" which are symmetric. We also assume the properties 

C{z)e : e > ? 7 |ep, C'{z)e : e > 0 (8) 

for all e G z G [0,1] and a constant p > Q independent of e and z. 

Furthermore, we choose the mobility m G C'(R x [0,1];R+) and suppose that the chemical 
energy density dl G C^(R) can be decomposed into 

'I'(c) = 'I'i(c) -I- 'I' 2 (c) for c G R, 

where d'l, ^'2 G C^(R) with convex and d'l > 0. 

In addition, we assume the following growth conditions: 

|vI/'(c)|<C(l + |cp*/2)^ 

|vI/'(c)|<C(|c| + l) 

for all c G R. Moreover, the mobility function should satisfy 

Cl < m{c, z) < C 2 

for all c G R, 2 ; G [0,1]. Here, Ci, (72 >0 denote constants independent of c and 2 , and 2* is the 
Sobolev critical exponent. 

The damage dependent potential / entering equation 631 is assumed to be a function of 
Ci([0,l];R+). 


(9a) 

(9b) 

( 10 ) 
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3 Main results 

3.1 Notion of weak solutions and existence results 

In what follows we define for fc > 1 the spaces 

:= {u e W^’P{n) I M > 0 a.e. in fl}, 

W'l’P{n) := {m € W^’P{n) I M < 0 a.e. in fl}, 

:= {u £ I M = 0 on Fd in the sense of traces}. 


Let the following initial-boundary data and volume forces be given: 

boundary data: b€ H\0,T] nW^’\0,T; 

initial values: c° £ H\n), £ 77^(17;R"), 

£ W^'P{Vl) with 0 < < 1 a.e. in t2, 

external volume forces: I £ L^(0, T; K")). 

A weak formulation of system (U)-(in]) is given in the following definition. 

Definition 3.1 (Weak solution) A weak solution of the PDE system (m)-® for the data 
(Z, 6, c°, z;°, z°) is a 5-tuple (c, u, z,/i, satisfying the following properties: 

• spaces: 

c £ L°°(0, T; H\n)) n iF^(0, T; {H\n))*), 

u £ L“(0, T; R")) n T; L^{Vl- R”)) n H^{Q, T; (fl; R"))*) 

with u = b on Fd x (0,T), u(0) = a.e. in O, dtu{0) = a.e. in ft, 

z £ L°°(0, T; W^’P{n)) n h\o, t- L^{n)) 

with z(0) = in z > 0 a.e. in fix, dtz < 0 a.e. in fix, 
p€L\0,T-H\n)), 
f £ L^{0,T;L\n)). 

• for all C £ L'^{0,T;H^{n))nH^i0,T;L^{n)) with ((T) = 0.- 

f {c — c°)dtC dx dt = f m{c, z)\7 fi ■ V( dx dt (11) 

J fix' fiT 

• for all C £ L‘^{0,T;H^{fl)) and for a.e. t £ (0,T): 

f pfdx= / (Vc • VC-t-Wc(c, e(M), z)C +'E'^(c)C) dx (12) 

Jn Jn 


for all C £ Ltrn(il;K") and for a.e. t £ (0,T); 


{dttU,C)m+ [ We(c, e(u),z) : e(C)dx = f l-fdx 

J Q J Q 


(13) 
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• for all C £ W^’^{n) and for a.e. t £ (0,T); 

0 < [ (iVz^-^Vz • VC + (VF,,(c, e{u),z) + f{z) + dtz + OC) dx (14) 

Jn 

• for all C £ and for a.e. t £ {0,T): 

0 > /" C(C ~ -2) dx (15) 

Jn 


• total energy inequality for a.e. t £ (0,T); 

£{t) + IC{t) + V{0, t) < £(0) + /C(0) + >Vext(0, t) 

with 


(16) 


free energy: Eft) J 

+ f {fiz{t)) + 'it{c{t)))dx, 
Jn 


-\yz{t)\P + i|Vc(l)p + W{c{t),e{u{t)),z{t)) ) dx 

p 2 ’ 


kinetic energy: 

^(i) := ^|9tu(t)pdx, 

dissipation: 

'D{0,t) := ( ( + m(c, z) V/rp) dxds, 

J n* 

external work: 

>Vext(0,t):= [ W^e{c,e{u),z) : e{dtb)dxds 


— / dtu ■ dttb dx ds / I ■ {dtu — dth) dx ds 

J Oi J Ot 

— f ■ dtb^ dx + f dfuit) ■ dibit) dx. 

J n Jo. 


Remark 3.2 Let z^ yL,ff) be a weak solution. Furthermore, if additionally 

c£ H\d,T-H\n)), u£ H^{0,T-H^{n-MF)), z £ H^{0,T-W^’P{n)), 
then for a.e. t £ (0,T) 

Zt - ApZ + W^zic, e(u), z) + f{z) + C + </> = 0 in {W^’P{n))* , 

C £ ^^wf^in) (^)> 

F £ 9dy^i,Pf^Q-j{dtz). 

Moreover, the energy inequality (HU) becomes an energy balance. 

The main aim of this work is to prove existence of weak solutions in the sense above. 

Theorem 3.3 Let the assumptions in Section\^be satisfied. To the given data I, b, c^, u^, v^, 
z^, there exists a weak solution of system O-dH) in the sense of Definition \3.1\. 
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3.2 Notion of weak solutions for a regularized system and existence 
results 

We will first study a regularized version of our phase separation-damage model. The passage to 
the limit is performed in Section [4.21 The regularization is needed in the existence proof in the 
first instance to pass from the time-discrete to the time-continuous system. 

The regularized PDE system for (5 > 0 is given by 


Ct = div(m(c, z)V^), 

= -Ac -I- Wc(c, e(u), z) -f 'I''(c) -I- 5ct, 
utt - div (We(c, e{u),z)) -I- SAu = I, 
zt - Apz + lT,z(c, e{u),z) + f{z) -f ^ -I- tp = 0, 
C e 5 /[o,oo)(2), 

‘P G 91{-oo,0]{zt), 


where the linear operator A : —>■ R"))* is defined as 


{Au,v)h‘^ := [ (V(V-«), V(Vu))i 
Jq 


dx := 




f d^Mfc d^Ufc 
Iq dxidxj dxidxj 


dx. 


A weak formulation of the regularized system such as in Definition 13.11 can be obtained with the 
corresponding modifications including the (5-terms. 

Definition 3.4 (Weak solution of the regularized system) A weak solution of the regular¬ 
ized PDE system for the data {l,b,c^,u'^,v^,z^) is a 5-tuple {c,u, z, p,,f) satisfying the following 
properties: 


• spaces: 

c G L°°{0,T; H\D)) n H\0,T- L'^{D)), 
with c(0) = c*^ a.e. in D, 

u e L°°(0, T; R”)) n T; R”)) n H^{0, T; (fl; R"))*) 

with u = b on Td x (0,T), u(0) = a.e. in D, dtu{0) = a.e. in D, 
z e L°°{0, T; n H\0, T; L‘^{n)) 

with z(0) = z^ in D, z > 0 a.e. in Dt, dtZ < 0 a.e. in VLt, 
fxGL^{0,T-,H\n)), 

^ e L°°{ 0 ,T;L\n)). 


• for all f G and for a.e. tG {0,T): 



{dtc) f dx dt 


/ m(c, z)V/r • VC dx dt 

J Ot' 


(17) 


for all C G H^{n) and for a.e. tG (0,r).' 


/ ^Cdx = f (Vc • VC + W^c(c, e(u),z)C + 'k'(c)C + S (dtc) C) dx 

Jfl JQ, 


(18) 
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• for all C £ (ri; M") and for a.e. t G (0,T); 

(i9ttM,C)Hi+ [ W^e{c,e{u),z) :e{C)dx + S{Au,C)m = [ l-(dx 

JQ Jci 

• for all C £ and for a.e. t £ (0,T): 

0 < f {\Wz\P-^Vz-WC + {W,,{c,e{u),z) + f'{z) + dtz + OC) dx 

Jn 

• for all f £ L“(r2) and for a.e. t£ (0,r); 

0 > f — z)dx 
Jn 


( 19 ) 


( 20 ) 


in 

• total energy inequality for a.e. t £ (0,T); 

£:(t) + ]C{t) + p(o, i) < £:(o) + /c(o) + Wext(o, t) 

with 

free energy: £{t) ■.= J (^^\Wz{t)\P + ^\Vc{t)\'^ + W{c{t),e{u{t)), z{t))^ dx 


( 21 ) 


( 22 ) 



+ + 4'(c(t))) dx -b ^{AUrft), Ur{t))H-^ 

kinetic energy: 

K^{t) := i|5tii(t)pdx. 

dissipation: 

^(O,^) := / ( -b d i9tcp-b m(c, z) V/rp) dxds, 

J fit 

external work: 

Wext(0,t) := f W^e{c,e{u),z) : e{dtb)dxds 

J r2t 


+ S [ {Au{s),dtb{s))H^ ds 

— / dtu • dtth dx ds / I • {dtu — dth) dx ds 

J J 

— / • dtb^ dx+ [ dtult) ■ dtblt) dx. 

Jn Jn 


The proof of the main result, see Theorem 13.31 is based on the existence of weak solutions for 
the regularized system. 

Theorem 3.5 Let the assumptions in Section\^be satisfied. To the given data I, b, c^, u^, v^, 
z^, there exists a weak solution of the regularized system in the sense of Definition \ 3. ^ 


4 Proof of the existence theorems 

4.1 Existence proof for the regularized system 

For the existence proof of the regularized system, we will use a semi-implicit Euler scheme solved 
by a recursive minimization procedure. 
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Let r > 0 denote the discretization fineness and let Mr := [T/t\ be the number of discrete 
time points. We fix a fc £ I,..., Mr and define the functional Tr ■ x x 

W^’P{n) ^ R by 


T^{c,u,z) := 




+ ^U,u)h^ + - 


|Vcp + W[c,e{u),z) + f{z) + 'I'(c) — l{kT) ■ u 


z — z. 


k-l 


L2 


r 

y 


U — 2Ur ^ + Ur 


k-2 


1 

y 


r.k-1 


Vo 


s 

+ y 


„fc-i 


L2 


1 dx 

2 

L2 


where Vb = {C £ = 0}. Note that the inverse operator A’^ ^ : Vb — 

Uo := {C £ /n Cdx = 0} of the operator : Uq ^ Vq given by 


{yu,m{Cr ^,Zr ^)V-)i2 

is well defined. The space Vb is endowed with the scalar product 

{u,v)vo ■■= {V{A-^u),m{c^-^,z’^~^)V{A-^v))L 2 . 

We refer to [GarOO] for details. 

A minimizer of J^r the subspace 


cGH^{V,)\ ( (c — c°) dx = Oder 

Jn 

x^z&W^’P{M)\d<z<z’^-^^ (23) 

obtained by the direct method in the calculus of variations is denoted by {Cr,Ur, Zr)- More 
precisely, by a recursive minimization procedure starting from the initial values {c^,u^,z^) and 
u~^ := we obtain functions (c^, z^) for fc = 0,..., Air- The velocity field Vr is set 

to {Ur — Ur~^)fT and br and Ir are given by b{Tk) and /(rfc). 

Let Wr £ {lr,br,Cr,Ur,Vr, Zr, Hr}j we introduce the piecewise constant interpolations Wr, 
Wr and the linear interpolation Wr with respect to time as 

Wr{t) := Wr with k = [t/r] , 

Wr (t) := with k = [t/rl , 

Writ) := Pw'; + (1 - k = [t/rl , P = ^ ~ ~ 

r 

and the piecewise constant functions tr and as 


I X {m £ I u\r^ = &(rfc)|r„} 


tr := [t/r] T = min{fcr | /c £ No and kr > t}, 
tp := max{0, tr — r}. 

We would like to remark that, by definition. Writ) = Writr) for all t £ [0,T] and 

_ Oyk-l I ^,k~'^ 

d,Vrit) = +Ur 

'T^ 
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for t S [t/r]. 

Since the functions z^) are minimizers, we obtain the following necessary conditions 

(Euler-Lagrange equations) by direct methods in the calculus of variations, cf. |HK111 IHK121 
IHK13a| : 

Lemma 4.1 There exists a time-discrete weak solution in the following sense: 

• spaces: 

Cr, Cf e L“(0, T; H\n)), c, e T; H\n)), 

Ur,Vr G Ur,Vr G T; 

G L°^{0,T;W^'P{n)), Zr G W^'^{0,T;W^'P{n)), 

f,r€L°°{0,T;H\n)), 

with 

Ct-(O) = c° a.e. in 11, Ut(0) = a.e. in H, Zr{0) = z^ in H, Ut(0) = a.e. in H, 

Ur = br on Ed x (0,T), Zr > 0 a.e. in fix, dt'zr < 0 a.e. in fix, 


: dt. 


(24) 


. for all C G L‘^{0,T-H^n)): 

/ (i9tCT-)Cda;dt = — / to(c“, z“)V^t-■ VCda: 

J J 

• for all C G iJ^(n) and for a.e. t G (0,T); 

f fJ.rCdx= f {VCr-VC + W^c{Cr,£iUr),Zr)C-\-'^'{Cr)C-\-6{dtCr)C)dx, (25) 

JQ JQ 

• for all C G (ri; K") and for a.e. t G (0,T); 

f dtVr • c dx + f W^eicr,e{ur), Zr) ■ e(C) dx + S{Aur, ()= [ C' fdx, (26) 
JQ JQ JQ 

• for a.e. t G (0,r) and for all C, G W^’P{V.) with 0 < C + Zr{t) < zf{t): 

0 < [ (\VZr\P-^VZr ■ VC + {W^Cr, e{Ur), Zr) + f'{Zr) + dtZr)C) dx. (27) 

jQ, 

Lemma 4.2 (A priori estimates) There exists a constant C > 0 independent of S such that 


(^) 

(ii) 


(Hi) 

(iv) 


IIVcT-||Loo(o,T;L2(n;R")) < C, ||9tCr||L2(o_T;L2(n)) < C, 

ll'*^T||L°°(0,T;ff2(n;R")) < C, || Xr || (0,T;L2(n;R'*)) < C, 

l|wT||L“(0.T;i?2(n;R"))nWi"”(0.T;L2(O;Rn)) < C, 

lkT||L~(0,T;L2(n;R"))nffi(0,T;(if2^(n;R"))*) < C'; 

IIVzT-||ioo(o_T;LP(n;R")) < C, \\dtZr\\L^(o^x-L^(Q)) < C, 

II V^T||L 2 (o,T;L 2 (n;Rn)) < C, ||m(c“, 27 )^/^V/r-r||L 2 (o^T;L 2 (n;R")) < C. 


10 









Proof. We split the proof into two steps. We first prove the a priori estimates (i), (ii) and 
(iv) and then we deduce estimate (hi). 

First a priori estimates. Testing (1^^ with t/Tt, testing (pK|) with Cr — c~, testing (1261) with 
Ut — u~ — (bT — b~), and adding everything, yield 


Tiit) + T 2 it) + nit) + nit) + nit) < o 


with 


Ti(t) := f \/Crit)-Vicrit)-c-it))dx+ [ V(Vu,(t)), V(V(u.(t) - u"(t)))) ds; 

Jn Jn 

+ / dt%it)-iurit)-u~it))Ax, 

Jn 

T2(t) —rf micrit),ZTit))\VpTit)\^ dx + T f 5 |dtCr(t)|^ da;, 

Jn Jn 

nit) := [ W,c(cT(t),e(uT(t)),Zr(t))(cr(t) - c“(t))da; 

Jn 

+ / W,e(cT(t),e(u.r(t)),2-r(t)) : e(wr(t) - Mt (t))da; 

Jn 

Tiit) := r f d>'(cT(t))dt&r(t) da; — r ( Fit) ■ dtu-rit) dx 
Jn Jn 


dt&T(t) + W,e(cT(t),e(Mr(t)),2-r(t)) : e(dt&7 


nit) ■■= -r J (dtVTit) 

-T J (_^5{\/i\/Urit)),\/iVidtbrit))))^„Xr.xn - Fit) ■ dtbrit)^ dx. 


In 

These terms are estimated in the following. 

• Convexity estimates yield 

1„_ ...,,2 1, 


nit) > -||Vc.(t)||i.(^) - -||Vc;(t)||i.(f,) 
s s 

+ n ll^(^'*^T(0)ll_L2(Q.][{nxnxn) ” ■^||V(Vu.,. (t) ) 11 ^2 (Q.Rn X n X n) 


2 ' 

+ 2 ~ (0llL2(n;R")- 


• We obtain for small rj > 0: 


T 2 it) >v[ (l|VMr(s)||i 2 (n;Rn) +<5||9tc^(s)|li2(n)) ds. 

J t-r 

• By the convexity argument and by Zr < , we gain 

W^eicTit),eiurit)),Zrit)) : e(Mr(i) - y-Fit)) 

> WiCrit),eiUrit)),Zrit)) - W(Cr (t) , £(«“ (t)), (0) 

> Wicrit),eiurit)),Zrit)) - W(c:;:(t), e(u“(t)), z“(t)) 

+ WncTis),eiuFis)),Zris))dtCris)ds, 
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and conclude (r; > 0 is chosen as small as necessary) 


T3{t)> [ {W{cr{t),e(urit)),Zrit)) -W{c^{t),e{u^{t)),z^ (t))) dx 

Jn 

+ / W^,c(Cr(s), e(MT(s)), ^r(s))9tcv(s) dxds 

Jt~ JQ. 

+ / Vh.c(Cr(s),e('U“(s)),Zr(s))9tCr(s)dxds 

Jt~ J vl 

> [ {W{cr{t),eiurit)),Zrit))-Wic~{t),e{u~{t)),z~{t))) dx 

Jn 

-Crt [ ||Vh,c(Cr(s),e(Wr(s)),^T(s))|li 2 (n) ds 

Jt- 

-Cv [ l|M^,c(Cr(s),e(Wr (s))>^r(s))||i2(n)ds 
-vj^_ l|5tCr(s)||i2(n)ds 

> [ {W{Cr{t),eiurit)),Zrit))-W{c~(t),e{u~{t)),z~{t))) dx 

Jn 

— Cn [ (||cr(s)||i 2 (Q) + ||c^ (s)||| 2 (f 2 ) + ||e('Ur(s))||| 2 (Q)) ds 




-Cr, 


K(s))lli2(n)ds-?7 l|5tCr(s)||i2(o) ds. 


Convexity of 'hi combined with growth condition (I9bll and Young’s inequality show 
T4{t)> f 'i/i{cr{t))dx - f 'i/i{c~{t))dx 

JQ. J Q 


“ + ll^r(s)|li2(n.Rn)) ds 

~ [ ^||'h2(''T('S))llL2(f2) + ||u,-(s)||^2(Q.Rr.)^ ds 

J It 

> [ 'i>i{cr{t))dx - [ dii{c~ {t)) dx 

Jn Jn 

- (l| 5 iC^(s)lli 2 (n) + || 4 (s)|li 2 (n.R„)) ds 

~ J^_ (ll'^T(®)llL 2 (n) + lkT(s)||L2(Q.Rn)^ ds. 


• By using the discrete integration by parts formula, i.e., 



dtVr ■ diW dx ds 


Vr{t) ■ dtbrit) da: 


(t) • dtbr{t 


t) dx 
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r / ^ris) 

j t-r J Q 


dtbr{s) - dtbr{s - t) 


da; ds, 


(28) 


we obtain 

T^{t)> - I Vr{t) ■ dtbr{t)dx + I v~{t) ■ dtbrit - T)dx 




(s)lli2(n) + Crj 


dtbr{s) - dtbris - t) 


L^{Q) 


ds 


(^I|cr(s)|li 2 (n) +?7l|e(Mr(s))|li2(n;Rnxn)) ds 
/ C'^l|e(5t&r(s))||i2(n;Rnxn)ds 

Jt- 

j_ (^^ll^(^'*^T('S))|li2(f2.Rnxnxn) + C?) 11 V ( V9t 6r (s) ) 11 ^2 (f2. J 

[ (^IKT(s)||i2(n;R™) + C'j?||5t&r(s)||L2(Q.Rr.)^ dS. 

J 


Summing over the discrete time points tr = 0,t, ... ^ kr for an arbitrary but fixed chosen fc G N, 
we can apply Gronwall’s inequality and obtain the following boundedness properties: 


ds 


II VCT-||ioo(o,T;L2(n;R")) <C, 
ll^tCr||L2(g < C, 

||V(VMT-)||L“( 0 ,T;L 2 (n;R"Xnxr.)) < C, 
||e(WT)||L~( 0 ,T;L 2 (n;RnXn)) < C, 
lkr||L'x,(0,T;L2(n;R")) < C, 

II V/aT-||i2(0,T;L2(n;R")) <C, 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 


where C > 0 is independent of r. Combining estimates dSH)-® with Korn’s inequality, we 
obtain 

l|WT||L=°(0,T;ff2(O.Rn)) < C. 

Consequently, by noticing Vr = dtUr, 

l|Wr||L“(0.T;_f/2(n;R'*))nWi.~(0.T;L2(n;R")) < C. 

A comparison argument in (1261) also gives 

lk'r||L~(o,T;Z,2(n;R"))nffi(O.T;(ffgj^(n;R’‘))*) < 

Second a priori estimates. Testing (E7)) with z~ — Zr, yields 

Zrit)\P-^V Zrit) ■ V{Zr{t) - ^-(t)) dx + ^t|| (t) || ^2 (fj) 

<-t[ (W^z{c{T{t)),e{ur{t)),Zrit))dtZrit) + f'{zr{t))dtZrit)) da;. 

Jn 
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Now we apply a convexity estimate and get 

< T?7||atZ^(t)|||2(o) +tC'^( 1 + ||c^(t)||i4(n) + ||e(Mr(i))|||,4(o.Rnx,.)). 

We end up with 

IIVZi-l|Loo(o,T;LP(n;R")) < C, 
ll^t^r ||L2(o,T;L2(n)) < C, 

where C > 0 is independent of r. 

By applying Poincare’s inequality, standard weak and weakly-star compactness results to the 
above a priori estimates, we obtain the following convergence properties. 

Lemma 4.3 (Convergence properties) There exist functions 

c G L^{0,T;H\n))nH\0,T;L‘^{Tl)), 

u G L°“(0,r;i?2(Q.K’^))n k'‘)) n T; (ff; M"))*), 

and subsequences (omitting the subscript) such that for all r > 1 and s < 2*; 


Cr,C^ -J> C 

weakly-star in L°°(0, T; i/^(n)). 

(35) 


strongly in W (0,T\ L‘^{Tl))^ a.e. infix. 

(36) 

C7- —y C 

weakly-star in 

(37) 


L°“(0, T; n (0, T; L^(0)), 

(38) 

Ur,Uf —>• M 

weakly-star in L°°{0, T ; K”)), 

(39) 


strongly in (0,T; H^{fl;W^)), a.e. infix, 

(40) 

Ut ^ U 

weakly-star in 



L°“(0, T; M’")) n Wi’°“(0, T; 1.2(0; 

(41) 

Vr,Vf: -)> dtU 

weakly-star in L°°(0, T; L2(0; R")), 

(42) 

Vr —>■ dtU 

weakly-star in 

(43) 


L^{0, T- 1.2(0; R")) n lli(0, T; (fl; R"))*), 

(44) 

Zt, Z“ ->• Z 

weakly-star in L°°(0, T; W^’^(n)), 



strongly in L^{0,T; L''{fl',W^)), a.e. infix, 

(45) 

^ Z 

weakly-star in 

(46) 


l.“(0, T; W^'P{fl)) n 1/1(0, T; L^{fl)), 

(47) 


weakly in L'^{0,T] H^{fl)), 

(48) 

m{cf,z~)^yfj,r to(c,z)2 V/r 

weakly in 1.2(0^ T; 1.2(0; R")) 

(49) 


as T \ 0. 
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Strong convergence of a subsequence of {Vz,-} in LP(r2T;R") can be shown as in [HK13a) by 
a tricky approximation argument. 

Lemma 4.4 (cf. |HK13a] ') There exists a sequence {rfcjfcgN such that Zr^ z in LP{0, T ; kk^>P(ri)) 
as Tk \ 0. 


For a time discrete solution of the regularized system, we can prove the validity of an energy 
inequality of type (ITBl) except the additional discretization error terms e:J.,..., which will turn 
out to converge to 0 in a certain sense as r 0. 


Lemma 4.5 (Discrete energy inequality) Let a time-discrete weak solution be given as in 
Lemma EH Then the following energy estimate is satisfied for a.e. t G (0,T); 


Sr{t)+ICr{t)+Vr{0,t) 


(4(5) + 4(s) + 4(5) + 4(5)) ds 


<f,(0)+/C,(0)+W,,ext(0,t) 


(50) 


with the discrete energies 

Srit) := Q|Vz^(t)|P + ^|Vc.r(i)P + W{Cr(t),e{Ur(t)),Zr(t)) fizrft))^ dx 

f s 

+ / 'i>{Cr{t)) dx -G -{AUr{t),Ur{t))H^, 

J o 

^r(i) := ^\Vr{t)\'^ dx, 

■Dr{0,t) ■= [ / (|9t2rP + i5|9tCrp + m(c“,z“)V/ir • V//T-)dxds, 

Jo Jo, 

Wr,ext(0,t) := [ [ W^e{Cr,e{Ur),Zr) : e{dtbT)dxds 


IQ Jo 

-\-S {Aur{s),dtbr{s))H^ ds 

Jo 

+ / / Ir ■ (dtUr - dtbr) dxds - 

Jo Jo ^ ! Jo 


dtUr — dtbr ) dx ds — / ■ dthr{d) dx 

o 

t-T 

Vt{^) ■ dthr{t) dx — 


[ Vt ( 5 ) 

Jo 


dtbris) - dtbr{s - t) 


dx ds, 


and the error terms 


4(0 := 


4(0 := 


4 ( 0 := [ 

Jo 


W(Cr(t),e{u,^ {t)),z,^ (t)) -W{Cr{t),e{u,^ {t)),Zrit)) 

T 

+ [ W^^(Crit),e(Ur{t)),Zr(t))dtZr(t)dx, 

Jo 

^ W{cf{t),e(uf{t)),zf{t)) - W{cr{t),e{uf{t)),zf{t)) 

+ [ W^c{Cr{t),e{Urit)), Zr{t)) dtCr{t) dx, 

Jo 


dx 


dx 


4'(c^ (0) - ^(Cr(0) 


dx + / dl'{Cr{t)) dtCr{t) dx, 
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r{t):= f 

JQ. 


f{Zr {t)) - fiZrit)) 


dx+ f'{Zr{t))dtZr{t)dx. 


Proof. We compute by using convexity of W with respect to e: 

/ W^e{Cr,e{Ur),Zr) ■■ e{Ur - U~)dx 
JQ. 

> / {W{Cr,e{Ur),Zr) -W{c~,e{u~),z~)) dx 
JQ 

+ / {W{cr,e{u~),z~)-W{cr,e{u~),Zr)) dx 
JQ 

+ / {Wic~,e{u~),z~)-W{cr,e{u~),z~)) da:. 
JQ 


(51) 


We test (E51) with Ur — u.^ — {br — b.^), apply (ICT) , use further convexity arguments and end up 
with 


\ Ikr(t)|li 2 - ^ \\v., + ^{AUr{t),Ur{t))H^ - ^{Au., (t ))//2 


+ f {W{Crit),e{Ur{t)),Zr{t)) -w{c.^{t),e{u.^{t)),z.^it))) dx 
JQ 

- [ dtVr{t) ■ {hr{t) - h~(t)) dx 

JQ 

+ / {W{cr{t),e{u-{t)),z-{t))-W{Cr{t),e{u-{t)),Zr{t))) dx 
JQ 

+ / {W{c-{t),e{u-{t)),z-{t))-W{cr{t),e{u-{t)),z-{t))) dx 
JQ 

< [ lr{t) ■ {ur{t) - U~{t) - {hr{t) - 6“(t))) dx 

JQ 

+ [ W^e{cT{t),e{ur{t)),Zr{t)) : e{br{t) - b~{t)) dx 
JQ 


(52) 


+ 5{Aur{t),br{t) - b^ 
Using the convexity estimate 


(53) 


[ • V{Zr - Z-) dx > hvZrWl, - -||Vz-|| 

JQ P P 


-IIP 
LP 


and testing (l?7l) with — z,- yield 


-iiv^.(t)iii. - k^z-mi^+rWdtMml^ 

p p 


< [ {W^z{e{Urit)),Zr{t))+ f{Zr{t))){z^{t)- Zr{t))dx. (54) 

JQ 

Next we test equation (|24ll with and m with {ct—c.^ ) and add the two derived equations. 
We obtain by means of the convexity property 

^ Vc, • V(c, - c;)dx > i||Vc.||i. - i||Vc;||i. 
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the estimate 


^l|Vc^(t)||i2 - ^l|Vc^ (t)||i2 + J^(w,c{Cr{t),e{Ur{t)),Zr{t)){Cr{t) - {t)) + 

'i>' {Cr{t)){Cr{t) - C~{t)) +Tm{c~{t),Z~{t))\7fir{t) ' ^ dx + 5 t|| 5* CV (t) < 0. (55) 


Adding the estimates (1551) - (l55l) . we end up with 


\\Vr(t)\\^2 - - ||x^ (t)||^2 + ^{AUrit),Urit))H^ - 

+ ^iivc.(t)iii. - \\\vc-mi. + kyzrimi, - Vz-(t)||i, 

+ T(^\\dtZrit)\\]^2 + S\\dtCr{t)\\l2 + j TO(c“(t),z“(t))V^r(t) ■ V/^r(t)dx 

J Q 

- [ dtVrit) ■ {br{t) - h~{t)) Ax 

Jn 

+ f {W{Crit),e{Ur{t)),Zr{t)) -W{c~{t),e{u~{t)),z-{t))) dx 

Jn 

+ [ f{zT{t))dx- [ f{z~{t))dx+ [ d>{cr{t))dx 

Jn Jn jn 

- / «'(c;:(t))dx + T(e^(t) + e^(t) + e^(t) + e)^(t)) 

Jn 

< [ lr{t) ■ {Ur(t) - U~{t) - {brit) - b~ {t))) dx 

Jn 


+ [ W^e{cr{t),e{ur{t)),Zr{t)) : e(6^(t) - b^ {t)) dx 

Jn 

+ S {Aur (t) , br {t) - b~ 

with the error terms e^{t), e^(t), e^{t) and e^(t). Summing over the discrete time points and 
taking into account the discrete integration by parts formula (1281) . we finally obtain the claim. 


Proof of Theorem \S.5[ We are going to establish the equalities and inequalities of the weak 
formulation m-m- 

• (Cahn-Hilliard equation) 

Because of the convergence properties (EH), dsni), (USD and (05]) we may pass to the limit 
in dn and obtain (Ha. 

To establish (HU, we first integrate (l25]l over time from t = 0 to t = T. The growth 
condition (1^ and the convergence properties (H51) . (1551) . (IMl) . (HOI) . (H5|) and (157)) allow us 
to pass to the limit in the integrated version of (1551) which shows (1181) . 

• (Balance equation of forces) 

By using the canonical embedding T^(fl;R”) ^ (i/p^^(fl;R”))*, it follows for all f G 

[ dtVr{t) • (dx = {dtVr{t),C,)H-^- 

Jn 

Keeping this identity in mind, integrating (1261) from t = 0 to t = T and passing to the limit 
r \ 0 by using (1551) and (1551) . (HiH) . (1551) and (155|) . we obtain (1T51) . 
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• (Variational inequality for z) 

To obtain the variational inequalities o and (ED), we can proceed as in |HK13a] . In 
particular, ED is valid for the subgradient 


C = -X{z=o}max|o,lV,^(c,e(u),2:) + /'(2)|, 

which satisfies ED, where X{z=o} is the characteristic function of the set {z = 0}. 
(Energy inequality) 

To treat the energy inequality (IHII)) . we set 
Ar{t) := 

/ + i|Vc^(t)P + W{c{t), e{Ur(t)), Zrit)) + f{Zr{t)) + T(Cr 

Jq 2 

- ^ + IV(c°, e(u°),z°) + /(z°) + T(c°)) dx 

+ f ^\vT{t)fdx- j ]-\v°\^dx+^{AUr{t),Ur{t))H^ -^{Au°,vP)h'^ 

«/ O «/ o 

— f Vr{t) ■ dtbr{t) dx + f v°-dtbT{0)dx 

Jq Jq 

Br{t) := 

J J (^\dtZr\^ + S\dtCr\^ + m{Cr , Zr)V dx ds 

— J y ■ {pt^T ~ 9tbPj dxds 

r-tx 


(56) 


10 Jn 

rt^ 


W^eiCT,e{Ur),Zr) : e{dtbr)dxds - S / {AUris),dtbris))H2 ds 


+ 


Vr (s) • 


dtbris) - dtbris - t) 


dx ds, 


/o Jn 

l-t-r 


El{t):= / ei(s)ds, El{t) ■.= / e^(s)ds, 

Jo Jo 

El{t):= [ e^(s)ds, E^{t) := f e^(s)ds. 

Jo Jo 

Then, (ISOl) is equivalent to 

Arit) + Br{t) + E^{t) + E^(t) + Ept) + Ept) < 0. 
Furthermore, by the a priori estimates, we observe that 

\Ar{t)\ + \BAt)\ + \El{t)\ + \E^At)\ + EUt)\ + \Ept)\ < C 


dx 


(57) 

(58) 


for all t € [0, T] and for all r > 0 (along a subsequence r^). Next, we consider the liminfT-\,o 
of each term in (ED separately. 


— By the already proven convergence properties and by lower semi-continuity arguments, 
we obtain 

rt2 pt2 

liminf / Ar{t) dt> A{t) dt for all 0 < ti < t 2 < T, (59) 
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where A is defined as A^- but c,-, Ut, z^, Vr and 6,- are substituted by their continuous 
limits. Note that this lim inf-estimate does not necessarily hold pointwise a.e. in t 
because, for instance, we do not know Vr{t) —>■ v{t) weakly in L^(r2;]R") for a.e. t (see 
621 )). 

Let 0 < < ^2 < r be arbitrary. By Fatou’s lemma, by (HH)) and by a lower 

semi-continuity argument, we obtain 


rt2 ptT n 

lim inf / / / 

’"NO Jq Jq 


Analogously, 


|9tZi-(s)P dxdsdt > [ (^liminf [ f |9tz^r(s)pdxds^ dt 

Jti V Jq Jq J 

*2 /-t 


> 


f [ |9tz(s)pda:dsdt. 
</ 0 o 


(60) 


pt2 pt-r p 

lim inf / / / 

’"NO Jf^ Jq Jq 


(5|9tCr(s)P dxds dt > 


('^2 pi 


(5|9tc(s)p dccdsdt (61) 


' ti «/ 0 


and, by (H^ . 


pt2 pt-r p 

lim inf / / / 

Jq Jq 


m{c^ (s), { 3 ))^iIt{s) ■ V^ir{s) dec ds dt 

pt2 pt 


> 


n m{c{s), z(s)) V^(s) • V/r(s) dec ds dt. (62) 
! 


Jti Jo J Cl 

Taking also (15811 and the already known convergence properties into account, we obtain 

Ct2 /'t2 


CC2 I't2 

liminf / Br{t)dt> / B{t)dt, 

’■NO 7*1 7*1 


(63) 


where B is defined as Br but Cr^Cr, Ut, Ut, , Zr^ 'zr, Jir and br are substituted by 
their continuous counterparts and ^tNd)-gtbx(t-T) dttb{t). 

Due to the differentiability of C we have 

C{z~) = C{zr) + C'{zr)iz~ - Zr) + r{z~ - Zr), -)■ 0 as *7 -)> 0. (64) 


Hence, we obtain 

/■* r 1 C(z-) - CjZr) 

7o Ja 2 T 


10 J{Zr {s)^Zr{s)} 

Because of 


[e{Ur) — e*{c)) : [e{Ur) — e*{c)) dxds 

'Av,a.,--.,)..--.. 


Zt Zt 


r(Zr - Zt) 

< 

M 

1 

1 

M 

L°°({z~=itZr}) 


: {e{Ur ) — e*(c)) dxds 
C{Zr)-C{Zr) 


(eK) - e*(c)) 
(65) 


Zr - Zr 


L’=°{{Zr ^Zr}) 

C'(zr)^^^ 
Zr - Zt 


<c. 
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and 


r{z^ -z-r) 


I^T -Zt\ 

gence theorem 


0 a.e. in as r \ 0 we conclude by Lebesgue’s generalized conver- 


r{z^ - Zr) 


Zr - Zr 


0 for every q > 1. 


L^iiZr ¥^Zr]) 


Using this and the already known convergence properties, we end up with 

left hand side of (EH) —t / W^z{c,e{u),z)dtzdLX<ls 

Jot 

and, consequently, —>■ 0 as t \ 0. Together with the uniform boundedness (IMl) . 

this implies 


rt2 

/ (t) dt —?► 0 as T \ 0 for all 0 < ti < t 2 < T. 

Ju 


The convergence 


pt2 

/ E^{t) ds —>■ 0 as r \ 0 for all 0 < ti < t 2 < T, 

Jti 


( 66 ) 


(67) 


can be shown as above. 

Noticing the linearity of e*, a short calculation yields 

/■*" f W{cf,e{uf),zf)-W{cr,e(uf),zf) 


dx 


/o Jn 


rt-r 


)(eK) - ■ e*idtCr)dxds. 


JQ Jn 

Due to the already known convergence properties, we obtain 
/■*" f W{c-,e{u-),z-) -W{c,e{u-),z-) 

Jo Jn ''' 


dxds—>■— / W^c{c,e{u),z)dtcdxds 


t Qt 


and, consequently, E^{t) —>■ 0 as r \ 0. Together with the uniform boundedness 
this implies 

pt2 

/ (t) dt —0 as T \ 0 for all 0 < ti < t 2 < T. 

Ju 


The claim 


pt2 

liminf / E^(t) dt > 0 for a\\ 0 < ti < t 2 < T 


( 68 ) 


(69) 


can be shown by the following arguments: On the one hand, convexity of Ti yields 

^l(Cr)-^l(Cr) ^ > 0. 

T 

On the other hand, by using the differentiability property of ^> 2 , we obtain (cf. (l64ll l 

„i,h IM ^ 0 as, ^ 0. 
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In the non-trivial case — c-r ^ 0, we can argue as follows. Since = 

’•(Ct -Ct) -Ct since is bounded in L^frir), it remains to show 

Ct -Ct t t V j yj 

^ —>• 0 in L^{Qt) as r \ 0. (70) 

Cr Cf 

Indeed, it converges pointwise to 0 a.e. in and applying the mean value theorem 
yields (here ^ G [min{c“,Cr},max{c“,Cr}]) 


r{Cr - Cr) 


42 (c-) - 42 (c^) 

Cr Cr 


Cr Cr 


< 

'1^2(01 + I^2(C.)| 


< C{1 + 1^1 + |ci-|) 

< C(1 + jc.,. I + 2 |ct|). 

Therefore, the left hand side is bounded in L“(0,T;L^ (^))- Lebesgue’s generalized 
convergence theorem yields dZOl). We end up with 

liminfT-\^o > 0 as t \ 0. Fatou’s lemma shows the claim. 

If we combine (15^1) . (1551) . ((551) . (1551) . (1571) and (1551) . we finally obtain 

0 > liminf f {Ar{t) + Br{t) + i?^(t) + E^{t) + E^{t) + E^{t)) dt 

> f {A(t) + B{t)) dt. 

Jti 

for all 0 < <1 < ^2 < F. Thus, A{t) + B{t) < 0 for a.e. t € (0, T) which is the desired energy 
inequality da. 

Hence, we obtain existence of weak solutions in the sense of Definition 13.41 □ 

4.2 Existence proof for the limit system 

We now study the limit (5 \ 0. For each 5 > 0, we obtain a weak solution (cs,us, zs, ns,is) in 
the sense of Definition 13.41 

Lemma 4.6 (A priori estimates) There exists a constant C > 0 independent of S such that 

(i) ||c5||Loo(o,T;ffi(n)) < C, VS\\dtCs\\ L^(o^T-,L^(n)) < C, 
l|5tC5||L2(o_T;(//i(n))*) < C, 

(ii) lk<5||L“(0,T;i/i(n;R"))rUVi'Oo(o,T;L2(a;R")) < C, ■\/^||M<5|U~(0.T;i/2(n;R")) < C, 

ll^<5||if2(0_T;(ff2j^(a;R"))*) < C, 

(Hi) |h5||L°°(0,T;VKi'P(n))nHi(0,T;L2(a)) < C, 

(iv) ||/i5||L2(0,T;ifi(a)) < C, ||m(c5, 25 )^/^V/i 5 ||i 2 ( 0 ,T;L 2 (O;Rn)) < C. 
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Proof. From the energy inequality (0^ . we infer the second inequality of (i), the first two 
inequalities of (m), (in) and the second inequality of (iv). By considering (HU), we get 

(dttusit), C)i/2 < c(\\e(u5(t))\\L2 + ||c5(t)||i2) l|e(C)llL= 

+ d||v (VM5(t)) IU 2 IIV (VC) IU 2 + 11 /IU.IICIU 2 

and, therefore, 

l|w5||ff2(o,T;(i/2^(n;R"))*) < C. 

Due to f^cs(t)dx = const, and the boundedness of ||VQ(t)||i 2 (f 2 ), we derive by Poincare’s 
inequality the first inequality of (i). 

From (dzj and (ED we obtain boundedness of fis{t) dx. Since || V/i 5 (f)||i 2 (f 2 ,j,) is also 
bounded, Poincare’s inequality yields the first inequality of (iv). 

Finally, we know from the boundedness of {V^^} in that {dtcs} is also bounded 

in L^(0,T; (iJ^(O))*) with respect to 5 by applying equation (fT71) . Hence, the third inequality 
of (i) is satisfied. 

Lemma 4.7 (Convergence properties) There exist functions 

ce L°°(0,T-,H\n))r\H\0,T;(H\n))*), 

u € L°°(0, T; H^(n- K")) n W^^^(0, T; L^(n; R”)) n H^(0, T; »”))*), 

z € L°^(0,T;W^'P(n)) n H^(Q,T; L^(n)), 

f,€L^(0,T-,H\n)) 

and subsequences (omitting the subscript) such that for all r >1 and s <2*: 


cs ^ c 

weakly-star in 



L^(0, T- H\n)) n H\0, T; (H\n))*), 

(71) 


strongly in LP'(TIt)^ a.e. in TIt, 

(72) 

Us ^ u 

weakly-star in L°°(0,T; H^(Tl]MT)), 

(73) 


weakly-star in VF^’°°(0, T; L^(H; R")), 

(74) 

zs ^ z 

weakly-star in 

(75) 


L°°(o,T; w^^p(n)) n H\0,T- L'^(n)), 

(76) 


strongly in L’'(0,T; L’'(H;R")), 

(77) 


a.e. in D,t, 

(78) 


strongly in LP(0,T;W^’P(n-R")), 

(79) 


strongly in C(Qt), 

(80) 

h-5 ^ h 

weakly in L^(0, T; iJ^(f2)), 

(81) 

—>■ m(c, z)^^^Vp, 

weakly in L^(0, T; L^(f2, R")) 

(82) 


as S \ 0. 

Proof. Lemma 321 reveals the existence of functions 
c€ L°°(0,T;ff\n)) 
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u G L°°(0, T; H^in- R”)) n T; L^{VL-, R")) n T; R”))*), 


z G L°°{Q,T-W^'^{n))f\H^{Q,T]L^{n)), 

^Ji&L\Q,T■H\^)), 
m(c,z)i/2VAi G L^(0,T; L2(f],R")) 

and subsequences indexed by 5k such that 

cs^ ^ c weakly-star in L°°(0, T; (83) 

usf. —>■ u weakly-star in (84) 

L°°(0, T; R")) n T; L'^{n-, R”")), (85) 

zsf^ ^ z weakly-star in L°“(0,T; Vl^^’^(r2)) n i/^(0, T; L^(ri)), (86) 

weakly in L^(0,T; i7^(0)), (87) 

m{cs^,zs^y^'^yfj, 5 k^w weakly in L^(0,T;L^(O;R")). (88) 


Due to the strong convergence properties of {c^^}, {zs^.} and the growth assumptions on the 
mobility function m, we infer 

w = m(c, 

In the following, we omit the subscript k. Furthermore, property (i) of Lemma 14.61 shows that 
{c^} converges strongly to an element c in L‘^{Qt) as i 5 \ 0 for a subsequence by a compactness 
result due to Aubin and Lions l |Sim86| L By choosing a further subsequence we also obtain 
pointwise almost everywhere convergence. 

By applying the same technique as for Lemmastrong convergence of Wzs in R”) 

can be obtained. Note that we need the assumption C'(-) > 0, see dHl). We conclude that 

ZS^Z strongly in LP(0,T;Wi’P(fI)). 

Furthermore, by Lemma 14.61 (iii), we find 

zs ^ z strongly in C(Dt) 

for a subsequence by an Aubin-Lions type compactness result (cf. [Sim86) L 

Next, we will proof our main result. 

Proof of Theorem \3.S\ . 

• (Cahn-Hilliard equation) 

Writing CZl) in the form 

/ (c5 — c°)i9tCdxdt = / m{cs,zs)V^JLs-VfdxAt, 

J J 

by only allowing test-functions C S L^(0, T; 7 L^(fI)) n 77^(0, T; L^(D)) with C,{T) = 0 we 
may pass to the limit by means of the convergence properties dZH), CID and m and 
receive (ED. 

Equation (ED can be obtained by integrating (1181) over time and taking advantage of the 
convergence properties (ED) (ED) (ED) (ED) (ED Lemmafib 
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• (Balance equation of forces) 

Integrating (fm from 0 to T and using (17^ , (ITHll . (IHHl) and the convergence S{Aus, 
—>■ 0 due to Lemma 14.61 (ii) we conclude 


(9ttu, C )//2 dt + / W^e{c,e{u),z) : e{C)dxdt = / l-(dxdt 


(89) 


for all C G L°°(0,T;iLpj^(n;]R”)). Therefore, (fOl) is true for all S 
K”) and a.e. t £ (0,T). Using the density of the set iLp^(fl;K”) in (here we 

need the assumption that the boundary parts Td and Tn have finitely many path-connected 
components, see |Berll| L we can identify dttu{t) G (iLp^^ (U; R”))* and (fUTl) is true for all 
C e i7jt^(U;R") and a.e. t G (0,r). Furthermore, dttu G 1^(0,T-, (iL^^(U;R"))*). 

• (Variational inequality for z) 

The variational inequality can be shown as in [HK13aj . We choose the following cluster 
points with respect to a subsequence: 


Xs ■= X{zi>o} X 

weakly-star in L°“(Ot), 

(90) 

Vs X{zs=0}n{W,^(cs,e(us),zs)+f'{zs)<0} V 

weakly-star in L°“(Ot), 

(91) 

Fs ■—X{zs> 0 }\l 2 e{cs))^F 

weakly in ^^(Ot’; R”^"), 

(92) 

5 •— X{zs=0}n{W^z{cs,e{us),zs)+f'izs)<0}^ 


(93) 

x^^^ieius)-e%cs))^G 

weakly in L^(Or; R"^"). 

(94) 


Note that since C'(-) is symmetric and positive definite matrix, its square root exists. By 
dza and (I5(I1) . we obtain for a.e. a; G {z > 0} 

X(a:) = 1, r]{x) = 0, F{x) = \l^^^^{e{u){x) - e*{cix))), Gix) = 0 (95) 

because of the following arguments: 

Let C S L^(f2T; R”^") with supp(C) Q {z > 0}. Then, by (IHHll . we obtain supp(C) C {zs > 
0} for all sufficiently small 5 > 0. By (19211 . we find 


: C dx dt —?► / F : C dx dt. 
J Qx” 


On the other hand, by (1751) . (note that <5 can be chosen arbitrarily small) 
/ F 5 :Cdxdt= [ \I {e{us) - e*(cs)) : (dxdt 


C'(z) 


(e(u) — e*(c)) : ( dx dt 


Thus, 


f \ — ^-^{e(u) — e* {c)) : ( dx dt = f F:(^dxdt 
Jn-r * 2 Jq 


The other identities in (IMl) follow analogously. 
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Now let C S L°°(0,T; W^l’^(O)). Taking (155)) into account, inequality (IT^ becomes by 
integration over time 


0 < / zs ■ y C + dtZsC) dt 

J Qt’ 


la 

+ 

+ 


/ 


{Z 5 > 0 } 


{W^z{cs,€{u5),zs) + f'{zs))Cdxdt 


{W^z(cs,eius),zs) + f'{zs))Cdxdt. 


J {zs=0}n{W^z{cs,e{us),zs)+f'{zs)<0} 

Applying lim sup, 5 x ^9 on both sides and multiplying by —1 yield 

0> lim [ (|V 25 |P“^Vz 5 • V(-C) + 9 t 2 ; 5 (-C)) dxdt 

+ liminf / (T' 5 )^(—dxdt + liminf f xs f {zs){—C) dx dt 
'5N0 Jot 

+ liminf j (G' 5 )^(—C) dxdt + liminf f r]s f'{zs){—()dxdt. 

Weakly lower semicontinuous arguments, the uniformly convergence property (IMl) and the 
properties listed in (l951l give 

0 > / {\yz\P-^yz ■ V(-C) + dtz{-0) dx dt 
J O-f 

+ [ {W,z{c,e{u),z) +f'{z)){-C)dxdt 

J{z>0} 

+ [ + G^) + {X + r])f'{z)) (-C) dx dt. 

J{z=0) 

This inequality may also be written in the following form; 

o< [ {\yz\p-^yz-yc + {w,zic,e{u),z) + f'iz) + dtz)c) dxdt 

+ [ {{F'^ + G'^) + {x + r])f{z)-W^z{c,e{u),z)-f{z))(:;dxdt. 

J |2; = 0]- 


Therefore, 


o< f {\yz\p-^yz-yc + iw,zic,eiu),z) +f{z) + dtZ + ^)C) dxdt 

with 

C := X{z=o}mm|o, {F^ + G^) + (x + ?? - l)/'(z) - kb,^(c, £(«), z)|. 
This proves da and da- 


• (Energy inequality) 

To prove the energy inequality da, we can proceed as in the proof of Theorem 13.51 
Integrating (1^ with respect to time on [ti,t 2 ] yields (0 < ti < t 2 < T) 



{-^s{t) + Bs{t) + Cs{t)) dt < 0 


(96) 
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with 


As{t) := 

J + ^|Vc5(t)P +W{cs,e{us(t)),zs(t)) + f{zs(t)) + ^'(c5(t))) dx 

-/ + + +f{z°) + 'i>{c°)^ dx 

+ [ ^\dtus{t)\^ dx - [ hv^l'^dx- [ dtus{t) ■ dtb{t) dx + [ v° ■ dtb° dx, 

Ja ^ Jn ^ Jn Ja 

Bs{t) := / [\dtzs\^+ 5\dtC5\^+m{cs,zs)^iJLs-'^lJLs)dxds 

Jo-t 

- / W^e{cs,e{us),zs) ■■ e{dtb)dxds + / dtus ■ dub dx ds 

J Q-t 

— l ■ {dtUs — dtb) dxds, 

Jo.t 

Cs{t) ■■= ^{Aus{t),us{t))H^ - ^{Au°,u°)h'^ - S J {Ausit), dibit)) h-^ dt. 

Let A be the corresponding integral expression to As, where cs, us and zs are replaced by 
c, u and z, respectively. Furthermore, let 

B{t) := / (jc^tzp + m(c, z)V/r • V^) dxds — / LF_e(c, e(u), z) : e(9tfe) dx ds 

J Qj J 

+ / dtu ■ dub dx ds — / I ■ {dtu — dtb) dx ds. 

J Ot J Ot 

The limit passage in (IMl) can be performed as follows. 

— Weakly lower semi-continuity arguments show 


pt2 pt2 

liminf / As{t)dt> / A{t)dt. 


Fatou’s lemma and weakly lower semicontinuous arguments for as well as the 
convergence property for cs, us, zs (see (ITTll . (17^ . (Tfl?!) . (175]) . (1771) 1 show (cf. (IHHI) - 


pt2 pt2 

liminf / Bs{t)dt> / Blt)dt. 
'5N-0 Jt^ Jt^ 


We have 


Cslt) > --{Au ,u )h^ —'5||M5(t)||//2(Q.Rn)||9t6(t)||jY2(n.Rn). 


By Lemma 14.61 liil. we obtain 


liminf [ Cs(t)dt>0. 

J. 


We end up with Alt) + Bit) dt < 0 for all 0 < ti < ^2 < IL- This proves (ITCl) . 
Putting all steps together, Theorem 13.31 is proven. 


□ 
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